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1.  Introduction 


It  has  long  been  known  that  chemistry  and  topology  are  both  important  in  the  glass 
forming  ability  of  metallic  melt.  The  primacy  of  these  features  are  reflected  in  2  of  the  3 
phenomenological  rules  1  for  formation  of  bulk  metallic  glasses  (BMGs)-  a  large  atomic  size 
mismatch  (topology)  and  a  large  negative  enthalpy  of  mixing  (chemistry).  While  these  rules 
are  rather  general,  recent  results  now  provide  specific  information  regarding  the  way  in  which 
topology  establishes  structure.  The  efficient  cluster  packing  (ECP)  model  ’  establishes 
representative  local  and  medium-range  structure  building  from  the  principle  of  efficient 
atomic  packing.  Topology  is  given  by  the  particular  sizes  and  concentrations  of  constituent 
atoms  in  the  metallic  glass. 

While  the  ECP  model  describes  how  topology  produces  structure,  it  cannot  yet  predict 
stability,  in  part  because  explicit  chemical  contributions  are  lacking.  If  structure  alone 
influences  stability,  then  topologically  identical  structures  might  be  expected  to  have  similar 
stabilities.  However,  it  is  well-known  that  metallic  glass  stability  in  a  given  alloy  can  be 
dramatically  altered  through  atom-for-atom  substitution  of  same-sized  atoms. 

The  complementary  roles  of  topology  and  chemistry  on  glass  forming  ability  can  be 
studied  by  careful  selection  of  metallic  glasses  that  are  topologically  identical  but  chemically 
distinct.  The  purpose  of  this  work  is  to  conduct  such  a  study.  Two  related  series  of  glasses 
will  be  studied;  the  topological  binary  Zrx(Cu,Ni)Y  series  and  the  topological  ternary 
ZrAAlB(Cu,Ni)c  series  of  glasses.  The  Zr-Al-(Ni,Cu)  system  is  chosen  for  study,  as  it 
represents  the  most  widely  studied  BMG  system.  Zr-Ni  and  Zr-Cu  both  form  binary  glasses 
with  solute  concentrations  in  the  range  of  30-50%  4.  While  Zr-Al  does  not  form  a  binary 
glass,  small  A1  additions  to  Zr-Ni  or  Zr-Cu  dramatically  improve  stability  5.  The  influence  of 
chemistry  (A1  substitution)  on  stability  is  analyzed  in  terms  of  chemical  bonding  and 
electronic  structure. 

A  structural  model  is  a  prerequisite  for  a  first-principles  study  of  electronic  structure 
and  chemical  bonding.  Either  experimental  or  calculated  (molecular  dynamics)  structures  can 
be  used  as  input;  here,  we  chose  to  start  with  the  experimentally  determined  crystal  structures 
for  CuZr2  and  NiZr2,  that  can  be  considered  as  parent  structures  of  the  glasses.  These 
compounds  have  the  appropriate  stoichiometry  for  extrapolation  to  some  glasses  under  study, 
and  their  crystal  structure  allows  the  calculations  for  extended  (periodic)  systems,  without 
artificial  boundary  effects.  The  Akmetal  substitution  can  be  easily  simulated,  and  a 
comparative  study  of  the  Zr-Cu-  and  Zr-Ni  systems,  with  or  without  substitution  by  Al,  was 
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achieved.  The  effects  of  the  Ahmetal  substitution  was  also  carefully  investigated  within  each 
system. 


2.  Structural  models 

The  experimental  crystal  structures  for  CuZr2  and  NiZr2  were  used  as  input  for  our 
calculations.  Al:Cu  and  Al:Ni  substitutions  in  CuZr2  and  NiZr?  were  simulated  by  using 
multiple  cells  (2x2x2  and  1x1x2,  respectively),  in  order  to  achieve  the  formulations  AlCu2Zr8 
and  AlNiiZi'8,  which  are  within  the  glass  forming  composition  range. 

In  CuZr2  with  tetragonal  structure  6  (space  group  I4/mmm,  2  formula  unit  (FU)  per 

o 

cell),  Cu  atoms  lie  in  a  Zr8  cubic  regular  site  (dcu-zr  =  2.848  A)  and  each  Zr  atom  has  four  Cu 

o 

first  neighbors  and  four  Zr  atoms  as  second  neighbors  (d/r./r  =  3.140  A).  The  crystal  structure 

n 

for  NiZr2  is  quite  different  (Al2Cu  structure  type,  space  group  I4/mcm,  4  formula  unit  per 

o 

cell);  Ni  atoms  lie  in  flattened  and  twinned  Zr8  cubic  site  (c/Ni-zr  =  2.761  A)  and  form  ID 

o 

chains  along  the  c  axis,  with  short  metal-metal  distances  (^Ni-Ni  =  2.634  A);  Zr  atoms  have 

o 

four  Ni  first  neighbors  and  three  Zr  neighbors  in  a  second  coordination  shell  (d/,./,  =  2.989  A 
xl,  3.079  Ax  2). 


3.  Computational  details 

Density  functional  theory  calculations  were  performed  using  the  augmented  plane 
wave  plus  local  orbitals  method  (APW+lo,  as  implemented  in  the  WIEN2k  code  "  )  and  the 
generalized  gradient  approximation  for  the  exchange-correlation  potential,  according  to  the 
formulation  by  Perdew,  Burke  and  Emzerhof  (PBE).11  Muffin-tin  radii  were  chosen  as  2.50 
a.u.  for  Zr  and  Cu,  and  2.49  a.u.  for  Ni.  1000  k-points  were  used  for  sampling  the  first 
Brillouin  zone.  In  the  case  of  2x2x2  and  1x1x2  multiple  cells,  100  and  500  k-points 
samplings  were  used,  respectively. 

Formation  enthalpies  were  calculated  by  subtracting  the  total  energies  for  elemental  A1 
(fee,  a  =  4.050  A),  Cu  (fee,  a  =  3.6144  A),  Ni  (fee,  a  =  3.5236  A,  ferromagnetic)  and  Zr  (hep, 

o  o 

a  =  3.233  A,  c  =  5.1475  A)  metals  from  the  compound  total  energy,  according  to  the 
stoichiometry.  Special  care  was  taken  for  calculating  all  energies  at  the  same  approximation 
level  (APW  +  lo  /  PBE),  and  with  identical  (or  closest  possible)  numerical  parameters 
(muffin-tin  radii,  cutoff  energies  ...).  In  all  calculations  for  elemental  metals,  both  spi- 
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polarized  and  non  spin-polarized  solutions  were  calculated,  to  ensure  that  the  calculated 
ground  state  has  the  correct  magnetization  (actually,  only  Ni  is  found  ferromagnetic). 

The  topologic  analysis  of  the  calculated  electron  density  was  performed  using  R.  F.  W. 
Bader’s  “Atoms  in  Molecules”  (AIM)  theory,12'15  as  implemented  in  the  WIEN2k  code.  The 
AIM  method  relies  on  a  non-ambiguous  definition  of  atomic  boundaries:  atomic  frontiers  are 
defined  as  surfaces  of  zero-flux  in  the  gradient  vectors  of  electron  density  (i.e.  ensemble  of 
density  local  minima  along  directions  passing  through  the  nuclei).  The  overall  electron 
density  is  thus  partitioned  in  well-defined  atomic  basins,  with  no  interstitial  space  left.  The 
atomic  charges  are  then  calculated  by  integration  of  the  electron  density  in  atomic  basins.  The 
AIM  method  constitutes  a  sophisticated  tool  for  the  characterization  of  chemical  bonding,  by 
the  analysis  of  features  at  critical  points  of  the  electron  density.  This  approach  will  be  detailed 
in  Part  4.3. 


4.  Results. 

4.1.  Formation  enthalpies 

Table  1  lists  the  energies  calculated  from  DFT  for  elemental  metals  and  the 
compounds  under  study.  The  formation  energies  deduced  for  CuZr2  and  NiZr2  are  in  very 
good  agreement  with  experimental  values  16,17  (-45.09  to  -49.78  kJ/mol  and  -95.24  kJ/mol, 
respectively).  This  is  a  key  point  for  checking  the  reliability  of  DFT  calculations  for  formation 
energies  in  these  compounds,  before  the  examination  of  substitution  effects. 

Al:Cu  and  Al:Ni  25%  substitution  in  CuZr2  and  NiZr2  induces  a  slight  destabilization 
in  the  case  of  CuZr2  (+1.3  kJ/mol,  i.e.  3%)  and  a  significantly  larger  one  in  NiZr2  (+16.7 
kJ/mol,  i.e.  17%).  This  might  been  due,  at  least  partly,  to  the  perturbation  of  strong  Ni-Ni 
interactions  by  the  Al:Ni  substitution  in  Ni  infinite  chains  along  the  c  axis. 
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Tablel.  Total  electronic  energies  per  formula  unit  (FU)  for  elemental  metals,  MZr2  (M  =  Cu, 
Ni)  and  their  Al-substituted  compounds,  and  formation  energies  calculated  from  DFT. 


Cu 

Ni 

Zr 

Al 

E  (Ryd/atom) 

-3310.05273 

-3041.65413* 

-7198.42269 

-485.64290 

CuZr2 

Al[/4Cu3/4Zr2 

NiZr2 

Ali/4Ni3/4Zr2 

E  (Ryd/FU) 

-17706.9293 

-17000.8259 

-17438.5743 

-16799.5588 

AHf  (Ryd/FU) 

-0.03119 

-0.03020 

-0.07479 

-0.06207 

Atf/(eV/FU) 

-0.42436 

-0.41086 

-1.01759 

-0.84454 

A/7/(kJ/mol) 

-40.94 

-39.64 

-98.20 

-81.49 

*Non-magnetic  Ni  metal:  E  =  -3041.64927  Ryd/atom. 


4.2  Electronic  structures 

From  our  DFT  calculations,  all  compounds  here  are  found  non-magnetic  metals.  The 
calculated  densities  of  states  (DOS)  with  contribution  of  atomic  muffin-tin  spheres  are  shown 
in  Fig.  1  and  2.  Figs  1(a)  and  2(a)  show  that  the  lower  part  of  the  valence  band  is  dominated 
by  Cu  and  Ni  3d-states,  with  broader  3d-bands  in  the  case  of  NiZr2,  due  to  strong  Ni-Ni 
interactions  (ID  Ni  ion  chains).  The  Zr  states  contribution  rises  with  energy  and  becomes 
dominant  at  Fermi  level.  Metal-Zr  interactions  occur  throughout  the  energy  range,  with  a  peak 
in  the  lower  part  of  NiZr2  DOS,  corresponding  to  bonding  Ni-Zr  interactions. 

The  Akmetal  substitution  induces  a  strong  perturbation  of  DOS  curves,  with  the  net 
stabilization  of  an  A1  3s-character  band,  with  bonding  metal- Al,  Zr-Al  interactions  (Figs  l(b, 
c)  and  2(b,  c)).  Thus  the  overall  chemical  bonding  is  strongly  affected  by  the  Akmetal 
substitution.  The  loss  in  formation  energies  with  the  substitution  is  the  resultant  of  stabilizing 
Al(3s)-metal  and  Al(3s)-Zr  bonding  interactions,  competing  with  antibonding  interactions  of 
same  character  (3  s  antibonding  bands  are  partly  occupied)  and  an  overall  destabilization  of 
the  electronic  structures  in  terms  of  translation  symmetry:  Al  ions  involve  a  different,  more 
ionic  chemical  bonding  than  Cu  and  Ni  ions,  breaking  thus  long-range  metallic  interaction 
chains. 
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Table  2  shows  the  atomic  charges  deduced  from  the  AIM  partitioning  of  the  calculated 
electron  densities.  NiZr2  appears  as  significantly  more  ionic  than  CuZr2.  The  net  negative 
charge  carried  by  A1  ions  is  only  63%  of  Cu  charge,  and  36%  of  Ni  charge,  acting  thus 
as  a  strong  perturbation  of  the  Madelung  potential  in  each  compound.  Ni  ions  bonded  to 
A1  ions  gain  a  significant  charge  and  lead  to  the  large  difference  in  A1  charge  between  the  Cu 
and  Ni  compounds.  Beyond  this  effect  of  short  Ni-Al  bonds,  the  differences  in  crystal 
structures  induce  a  distribution  in  Zr  charges  in  Ali/4Cu3/4Zr2,  while  a  single  Zr  site  (and 
charge)  is  found  in  Ali/4Ni3/4Zr2. 


Table  2.  Net  charges  calculated  by  the  AIM  partitioning  of  the  total  electron  density. 


Cu/Ni 

Zr 

A1 

CuZr? 

-1.08 

+0.55 

AI1/4CU3/4ZP, 

-1.08  /  -1.11 

+0.49  /  +0.50  /  +0.53 

-0.68 

NiZr2 

-1.31 

+0.65 

Alj/4Ni3/4Zr2 

-1.48/ -1.30 

+0.59 

-0.46 
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Figure  1.  Total  and  projected  DOS  for  CuZr2  and  Ali/4Cu3/4Zr2,  calculated  by  the  L/APW+lo 
DFT  method. 
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Figure  2.  Total  and  projected  DOS  for  NiZr2  and  Ali/4Ni3/4Zr2,  calculated  by  the  L/APW+lo 
DFT  method. 
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4.3  Chemical  bonding 


The  topology  of  the  electron  density  p(r)  is  characterized  by  critical  points,  where  the 
first  derivatives  of  p(r)  vanish.  Thus  critical  points  are  maxima  (nuclei),  minima  or  saddle 
points  in  p(r);  chemical  bonding  features  are  mainly  associated  to  the  distribution  and 
characteristics  of  saddle  critical  points.  A  critical  point  at  rc  can  be  further  characterized  by 
the  Hessian  eigenvalues  of  p  at  rc  (i.e.  laplacian  of  p(r)  along  local  principal  axes),  by  its 
signature  (algebraic  sum  of  the  signs  of  the  three  Hessian  eigenvalues),  by  the  ratio  0  =  A,// 
/  |AJ  of  the  laplacian  value  parallel  to  the  bond  path  to  its  main  perpendicular  component,  and 
by  the  value  of  p(rc)  itself.12'14 

The  nature  of  chemical  bonding  can  be  discussed  on  the  basis  of  critical  points 
characteristics: 14,16 

Shared  interactions  (covalent,  polar  bonding)  correspond  to  V2p(rc)  <  0,  0  <  1, 
large  p(rc), 

Closed-shell  interactions  (ionic,  hydrogen  bonding)  correspond  to  V"p(rc)  >  0,  0  > 
1,  relatively  low  p(rc), 

Metallic  bonding:  V2p(rc)  >  0,  low  p(rc), 

In  general,  cohesive  energies  vary  as  p(rc). 

These  limiting  cases  and  general  rules  are  illustrated  in  Table  3,  where  critical  points 
characteristics  for  diamond  (covalent  and  high  cohesive  energy)  and  sodium  chloride 
(archetypal  ionic  compound)  are  displayed. 


Table  3.  Critical  points  characteristics  for  diamond  and  sodium  chloride,  from  the  AIM 
analysis  of  the  calculated  electron  density.  All  values  are  in  atomic  units  (a.u.)  xlO4. 


Bond  path 

Signature 

Hessian  eigenvalues 

Trace 

p(rc) 

A,// 

>> 

1- 

V2p(rc) 

C-C 

-1 

1168 

-4211  -4211 

-7254 

2387 

Na-Cl 

-1 

752 

-89  -89 

575 

117 
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Table  4  lists  the  same  characteristics  for  critical  points  in  elemental  metals,  as 
calculated  in  our  study.  As  expected,  values  indicate  an  intermediate  situation  between 
covalent  and  metallic  bonding  for  Cu,  Ni  and  Al;  cohesive  energies  (p(rc))  vary  as  Ni  >  Cu  > 
A1  (in  agreement  with  tendencies  in  experimental  bulk  modulus),  correlated  to  a  partial 
transfer  from  covalent  to  metallic  bonding  in  the  series  (0  ,  V"p(rc)). 

The  case  of  zirconium  metal  is  more  complex,  with  large  0  values  and  moderate 
electron  density  at  critical  points,  a  situation  with  closed-shell  (ionic  cores)  and  metallic 
interactions.  The  main  type  of  critical  points  (Zr-Zr  1st  neighbors)  as  a  unusual  +1  signature, 
and  critical  points  are  found  for  2nd  and  3rd  neighbors  too,  indicating  a  complex  or  diffuse 
(metallic?)  chemical  bonding  picture. 


Table  4.  Critical  points  characteristics  for  elemental  metals,  from  the  AIM  analysis  of  the 
calculated  electron  density.  All  values  are  in  a.u.  xlO4. 


Bond  path 

Signature 

Hessian  eigenvalues 

Trace 

Z// 

Zx 

V2p(rc) 

Cu-Cu 

-1 

164 

-161 

-109 

-106 

391 

Ni-Ni 

-1 

108 

-199 

-136 

-226 

462 

Al-Al 

-1 

109 

-103 

-28 

-22 

299 

Zr-Zr 

+  1 

292 

35 

-52 

274 

257 

2nd  neigh. 

-1 

150 

-57 

-48 

44 

294 

3ld  neigh. 

+  1 

12 

12 

-7.5 

17 

202 

Some  tendencies  of  elemental  metals  are  retrieved  in  the  critical  points  of  the 
calculated  electron  density  for  CuZr2,  Ali/4Cu3/4Zr2  (Table  5),  NiZr2  and  Ali/4Ni3/4Zr2  (Table 
6).  Ni-Ni  interactions  are  stronger  than  Cu-Cu,  and  Zr-Zr  bonding  remains  complex  with  +1 
signatures  and  long-range  interactions.  All  traces  are  positive,  indicating  competing 
interactions  as  found  in  iono-covalent  and  metallic  compounds. 

Ni  ID  chains  in  NiZr2  contain  shared-interactions  critical  points,  similarly  to  Ni  metal. 
This  participates  to  the  large  stability  of  NiZr2,  despite  the  lower  dimensionality  of  Ni-Ni 
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interactions  than  in  the  metal.  The  overall  p(rc)  values  are  larger  in  NiZr2  than  in  CuZr2,  in 
agreement  with  the  larger  formation  energy  of  the  former. 

Cu-Zr  interactions  dominate  in  CuZr2  and  the  Cu-Cu  bonding  is  weaker  than  in  Cu 
metal;  the  crystal  structure  plays  an  important  role  here,  but  might  be  seen  itself  as  a 
consequence  of  weaker  Cu-Cu  bonding  compared  to  Ni-Ni.  V“p(rc)  and  p(rc)  values  in  Tables 
5  and  6  show  that  positively  charged  Zr  ions  (Table  2)  form  stronger  Zr-Zr  bonds  than  in  Zr 
elemental  metal. 

As  seen  in  Tables  5  and  6,  Al:Cu  and  Al:Ni  substitutions  induce  a  narrow  but 
significant  distribution  of  Hessian  eigenvalues  and  p(rc)  values,  indicating  a  perturbed 
chemical  bonding  distribution  in  the  two  compounds.  Critical  points  for  Zr-Zr  interactions 
appear,  that  were  not  present  in  the  parent  compounds,  showing  thus  a  higher  degree  of 
competition  between  interactions  in  the  Al-substituted  systems. 

Furthermore,  a  strong  perturbation  of  the  crystals  bonding  schemes  is  induced  by  the 
very  different  nature  of  Al-metal  bonds,  compared  to  the  parent  interactions  involving  Cu  and 
Ni.  This  appears  clearly  too  in  the  electron  density  maps  shown  in  Figures  3  and  4. 

A  point  that  remains  unclear  is  the  absence  of  Ni-Ni  and  Ni-Al  critical  points  for 
Ali/4Ni3/4Zr2  (Table  6),  whereas  the  strong  ID  interactions  found  in  NiZr2  still  exist.  This  is  a 
numerical  artifact  still  to  be  resolved. 

From  Tables  5  and  6,  the  perturbations  induced  by  the  A1  substitution  appear  more 
numerous  in  CuZr2  than  in  NiZr2;  however,  the  values  distribution  is  somewhat  broader  in 
Ali/4Ni3/4Zr2,  and  bond  paths  slightly  larger  than  bond  lengths  in  this  system  indicate  too  an 
instability  in  the  crystal  structure.  The  chemical  bonding,  critical  points  nature  and 
distribution  are  strongly  perturbed  and  disordered  in  both  compounds  by  the  Akmetal  partial 
substitution;  this  might  strongly  influence  the  tendency  to  glassy  transition,  together  with  the 
global  stability  of  the  parent  crystalline  systems  and  the  non  trivial  features  of  Zr-Zr  bonding. 
On  the  other  side,  Ni-Ni  bonding  appears  to  be  both  strong  and  directional,  even  in  binary  or 
ternary  systems;  this  may  be  an  unfavorable  factor  for  a  glassy  transition. 
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Table  5.  Critical  points  characteristics  for  CuZr2  (top)  and  Ali/4Cu3/4Zr2  (bottom),  from  the 
AIM  analysis  of  the  electron  density.  All  values  are  in  a.u.  xlO4  except  path  lengths  (A). 


Bond 

Path 

Signature 

Hessian  eigenvalues 

Trace 

p(rc) 

path 

length 

Z// 

K 

K 

V2p(rc) 

Cu-Zr 

2.85 

-1 

402 

-120 

-118 

165 

322 

Cu-Cu 

3.22 

-1 

271 

-42 

-40 

189 

246 

Zr-Zr 

3.14 

+1 

296 

63 

-67 

292 

281 

Cul-Zr 

2.85 

-1 

399 

-105 

-116 

178 

320 

-1 

401 

-105 

-116 

179 

320 

-1 

398 

-106 

-119 

173 

322 

-1 

399 

-107 

-116 

177 

321 

Cu2-Zr 

2.85 

-1 

404 

-117 

-121 

166 

323 

-1 

401 

-118 

-123 

160 

324 

-1 

403 

-117 

-120 

166 

322 

Al-Zr 

2.85 

-1 

276 

-80 

-86 

110 

332 

-1 

278 

-80 

i 

00 

115 

330 

-1 

280 

-82 

-83 

115 

331 

Cul-Cu2 

3.22 

-1 

270 

-42 

-40 

188 

246 

Al-Cul 

3.22 

-1 

207 

-48 

-46 

113 

275 

-1 

208 

-49 

-46 

113 

274 

Zr-Zr 

3.14 

+1 

300 

70 

-76 

294 

281 

+1 

297 

67 

-72 

292 

280 

+1 

302 

64 

-73 

293 

281 

+1 

296 

68 

-66 

298 

278 

+1 

304 

57 

-73 

288 

282 

Zr-Zr 

3.38 

+1 

156 

85 

-33 

207 

251 

+1 

153 

85 

-35 

203 

253 

Zr-Zr 

3.45 

+1 

128 

113 

-25 

215 

230 

+1 

130 

113 

-25 

218 

230 
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Table  6.  Critical  points  characteristics  for  NiZr2  (top)  and  Ali/4Ni3/4Zr2  (bottom),  from  the 

4  0 

AIM  analysis  of  the  electron  density.  All  values  are  in  a.u.  xlO  except  path  lengths  (A). 


Bond 

Path 

Signature 

Hessian  eigenvalues 

Trace 

P(rc) 

path 

length 

A,// 

V2p(rc) 

Ni-Zr 

2.76 

-1 

390 

-137 

■143 

109 

384 

Ni-Ni 

2.63 

-3 

-121 

-121 

-8 

-251 

430 

Zr-Zr 

2.99 

-1 

486 

-43 

-36 

407 

291 

Zr-Zr 

3.08 

+1 

398 

30 

-19 

409 

255 

3.56 

+1 

128 

25 

-16 

137 

225 

Nil-Zr 

2.76 

-1 

386 

-141 

■132 

112 

390 

2.77 

-1 

401 

-137 

133 

131 

378 

Ni2-Zr 

2.77 

-1 

397 

-144 

■137 

113 

382 

Al-Zr 

2.80 

-1 

274 

-125 

-85 

64 

387 

Zr-Zr 

2.99 

-1 

486 

-44 

-33 

410 

286 

Zr-Zr 

3.08 

+1 

398 

24 

-19 

402 

255 

3.13 

+1 

327 

70 

-31 

365 

264 

3.38 

+1 

197 

21 

-17 

201 

240 

3.38 

+1 

201 

9 

-21 

189 

242 

3.56 

+1 

121 

29 

-13 

137 

223 

3.75 

+1 

88 

37 

-13 

111 

234 

16 


Figure  3.  Electron  density  map  for  CuZr2  (top)  and  Ali/4Cu3/4Zr2  (bottom),  calculated  with  the 
L/APW+lo  DFT  method. 


Figure  4.  Electron  density  map  for  NiZr2  (top)  and  Ali/4Ni3/4Zr2  (bottom),  calculated  with  the 
F/APW+lo  DFT  method. 
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5.  Conclusion 


The  role  of  chemistry  on  the  enhancement  of  the  formation  of  glassy  phase  in  binary 
Cu-Zr  and  Ni-Zr  alloys  by  A1  substitution  was  carried  out  by  analysis  of  the  electronic 
structure,  chemical  bonding,  nature  and  distribution  of  critical  points  of  the  electron  density  in 
Cu-Zr,  Ni-Zr,  Al-Cu-Zr  and  Al-Ni-Zr  compounds  with  compositions  within  the  glass  forming 
range. 


It  was  found  that  the  overall  chemical  bonding  is  strongly  affected  by  Al:  metal  partial 
substitution  through  the  following  features: 

Al  ions  carry  a  much  lower  net  charge  than  Cu  and  Ni  ions,  and  break  the  site 
equivalencies  for  the  net  charges  of  Cu,  Ni  and,  in  the  case  of  NiZr2,  Zr.  Therefore, 
the  electrostatic  field  at  each  site  (Madelung  potential)  is  strongly  perturbed  and 
has  a  lower  symmetry. 

Al-metal  bonding  involves  Al  3s  orbitals,  that  have  a  different  symmetry  and 
extent  than  Cu  and  Ni  3d  orbitals.  Thus  Al  ions  break  metal  interaction  chains  and 
related  energy  dispersion  (bands)  and  destabilize  the  host  crystal  lattices.  Strong 
Ni-Ni  bonds  are  specifically  affected,  leading  to  a  significant  loss  in  formation 
energy  for  Al:NiZr2. 

The  distribution  of  critical  points  in  the  electron  density  is  strongly  affected  by  the 
Akmetal  substitution;  the  chemical  bonding  patterns  become  more  complex,  with 
an  increase  of  next-nearest  neighbors  interactions  for  Zr-Zr  bonds.  The  crystal 
structure  is  under  strain  due  to  competing  interactions  and  symmetry  lowering. 

These  results  suggest  that  the  enhancement  of  the  formation  of  the  glassy  phase  by  Al 
substitution  in  these  systems  has  a  strong  chemical  origin  and  can  be  attributed  to  the 
destabilization  and  perturbation  of  the  bonding  pattern  of  the  crystal  structure. 

Trends  found  in  this  study  should  be  transferable  to  other  crystal/glass  compositions  in 
this  family.  More  general  results  concern  the  strength  and  directionality  of  Ni-Ni  bonding,  and 
the  transformation  of  the  covalent  character  of  bonding  in  elemental  metals  into  a  much  more 
iono-covalent  character  in  binary  compounds.  The  agreement  between  calculated  and 
experimental  formation  energies,  and  trends  in  cohesive  energies  for  elemental  metals,  make 
this  approach  promising  and  reliable  for  further  stability  and  glass  formation  studies. 
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